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Abstract 

In  this  paper,  we  consider  the  problem  of  finding  locations  of 
several  new  facilities  in  an  imbedded  tree  network  with  respect  to 
existing  facilities  at  known  locations  so  as  to  satisfy  distance  con- 
straints, which  impose  upper  bounds  on  distances  between  pairs  of 
facilities.  It  is  known  that  the  existence  of  a feasible  solution 
to  the  distance  constraints  is  related  to  shortest  paths  through  an 
auxiliary  network,  which  has  as  arc  lengths  the  upper  bounds  on  pairwise 
facility  distances.  This  relationship  takes  the  form  of  necessary  and 
sufficient  conditions,  termed  the  separation  conditions.  We  relate 
"tight"  separation  conditions  to  the  solution  of  multifacility  minimax 
location  problems  and  efficient  solutions  to  multiobjective  multifacility 
location  problems.  In  addition,  we  provide  a proof  that  the  SLP  algorithm 
of  [8]  is  in  fact  an  algorithm  of  lowest  order  of  computation  for  deter- 
mining whether  or  not  the  distance  constraints  are  consistent. 
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1 . Introduction 

Network  location  problems  involve  locating  new  facilities  on  a 
network,  such  as  a transport  network,  with  respect  to  existing  facil- 
ities lying  on  the  network  at  known  locations.  For  such  problems  dis- 
tances between  facilities  are  of  interest:  these  distances  are  actual 
"shortest  route  distances"  as  determined  from  the  network.  Associated 
with  many  network  location  problems  is  a collection  of  distance  con- 
straints. These  distance  constraints  specify  known  upper  bounds  on 
distances  between  pairs  of  facilities  (pairs  of  new  facilities,  and  pairs 
of  new  and  existing  facilities)  so  that  the  facilities  will  not  be  too 
far  apart  from  one  another. 

Distance  constraints  appear  to  be  a recurring  feature  in  network 
location  problems.  In  the  problem  of  deciding  where  to  locate  a 
fire  station  [18],  it  may  be  relevant  to  place  restrictions  on  the 
location  via  distance  constraints  so  that  the  response  time  to  the  site 
of  any  potential  fire  may  not  be  too  large.  The  location  of  booster 
pumps  in  a pipeline  system  might  well  involve  distance  constraints  in 
order  to  locate  the  pumps  so  that  the  pressure  in  the  pipeline  does  not 
fall  below  a certain  level.  Also,  certain  minimax  location  problems  on 
networks  have  an  equivalent  formulation  in  terms  of  distance  constraints 
[8],  since  typical  objective  functions  in  such  problems  involve  dis- 
tances between  facilities,  permitting  equivalent  formulations  with 
distances  as  constraints. 

Our  work  extends  earlier  results  [8]  obtained  for  distance 
constraints  for  the  case  where  the  underlying  network  is  a tree,  that 
is,  it  is  connected,  undirected,  has  positive  arc  lengths  and  has  no 
cycles.  Such  networks  tend  to  occur  when  only  a "sparse"  network  is 


justifiable  because  of  limited  need  for  the  network  or,  alternatively, 
when  it  is  very  expensive  to  provide  redundant  links  (as  in  some  road 
networks).  In  addition,  some  pipeline,  and  associated  electrical  con- 
trol networks,  are  tree  networks  [3]. 

We  now  give  a brief  overview  of  our  main  results;  the  necessary 
notation  and  definitions  appear  in  Section  2,  including  the  "distance 
\ constraints."  Also  we  state  the  "separation  conditions",  necessary  and 

sufficient  conditions  for  a feasible  solution  to  exist  to  the  distance 
constraints.  It  is  shown  in  [8]  that  in  order  for  the  separation  condi- 
tions to  be  both  necessary  and  sufficient  for  the  existence  of  a feasible 
solution  to  the  distance  constraints  the  network  must  be  a tree.  Thus  in 
this  paper  we  restrict  ourselves  to  the  case  where  the  network  is  a 
tree.  In  Section  3 we  establish  preliminary  results  and  study  "tight 
paths",  paths  which  involve  an  auxilary  network  defined  by  the  distance 
constraints.  The  properties  of  Section  3 provide  the  basis  for  our 
subsequent  results . 

In  Section  4 we  give  necessary  and  sufficient  conditions  for  a 
feasible  solution  to  be  an  optimal  feasible  solution  to  a minimax 
location  problem,  and  present  an  example.  Minimax  location  problems 
occur  when,  for  example,  new  facilities  are  to  be  placed  so  as  to  mini- 
mize the  maximum  of  travel  times  between  given  pairs  of  new  facilities, 
and  given  pairs  of  new  and  existing  facilities.  Thus,  the  new  facilities 
are  to  be  placed  so  as  to  provide  quick  "service"  to  existing  facilities, 
while  at  the  same  time  being  able  to  provide  service  to  one  another.  We 
list  as  references  a number  of  papers  which  consider  minimax  location 
problems  on  networks  (See  references  [4],  [5],  [9],  [10],  [14],  [15]). 


J 
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In  Section  5 we  study  a multiobjective  location  problem  for  which 


f 
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it  is  desired  to  locate  certain  pairs  of  new  facilities  and  certain 
pairs  of  new  and  existing  facilities  as  close  together  as  possible. 
Efficient  location  vectors  to  such  problems  are  such  that  in  order  for 
any  new  facility  to  be  closer  to  some  facility  than  it  already  is,  it 
must  in  turn  be  placed  farther  from  some  other  facility.  This  parti- 
cular multiobjective  location  problem  is  similar  in  spirit  to  a problem 
in  the  plane  using  rectilinear  distances  studied  by  Wendell,  Hurter  and 
Lowe  [19],  and  Chalmet  and  Francis  [1].  Other  examples  of  multiobjective 
location  problems  on  networks  can  be  found  in  [11],  [12],  [13]  and  [16]. 

Also  in  Section  5 we  employ  the  idea  of  efficient  location  vectors 
to  demonstrate  that  any  algorithm  which  checks  the  distance  constraints 
for  consistency  (i.e.,  whether  or  not  the  distance  constraints  have  a 
feasible  solution)  is,  for  the  worst  case,  at  least  of  order  m + mn, 

where  m and  n are  the  number  of  new  and  existing  facilities  respectively. 

2 

We  further  point  out  that  the  "SLP"  algorithm  of  [8]  is  of  order  m + mn, 
so  that,  in  a well  defined  sense,  the  algorithm  is  in  fact  a "best" 
algorithm. 

Throughout  Sections  4 and  5 we  rely  heavily  on  properties  of  tight 
paths,  which  provide  the  underlying  unifying  principles  in  our  work. 


i 
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We  suppose  we  are  given  T,  an  imbedded  undirected  tree  having 
positive  arc  lengths,  such  that  the  distance  d(x,  y)  between  points 
x,  y £ T is  the  well-defined  shortest  path  distance  and  is  a metric, 


as  discussed  in  (5].  We  are  given  existing  facility  locations 

a^,  a^ aQ  in  T.  We  wish  to  find  new  facility  locations  x^,  x^ 

in  T,  where  x.^  is  the  location  of  new  facility  i,  i = 1,  m. 

As  in  [6],  given  points  x,.y  £ T,  we  define  the  line  L(x,  y)  as  the 
union  of  all  points  in  the  shortest  path  between  x and  y;  equivalently, 
L(x,y)  = {z  £ T:  d(x,z)  + d(z,y)  = d(x,y)}.  In  addition,  given  a col- 
lection of  points  P = {Pj,  P^}  contained  in  T,  we  let  H(P)  denote 

the  convex  hull  [6]  of  P,  the  smallest  (imbedded)  subtree  of  T 
containing  P.  Corresponding  to  the  n existing  facilities  and  m new 
facilities,  we  are  given  a set  of  Distance  Constraints,  which  for 
convenience  we  refer  to  as  DC,  as  follows: 


d(Xj , xk)  < bjk  , (j  , k)  £ Ig, 
dCx^  a^)  < c.j  , (i,  j)  £ Ic, 


(2.1) 


where  the  b..  and  c. . are  known  positive  constants.  We  remark  that 
it  need  not  be  the  case  that  the  set  Ig  includes  all  possible  pairs  of 
new  facility  indices,  nor  1^  includes  all  possible  pairs  of  new  and 
existing  facility  indices. 


Corresponding  to  DC,  we  define  Network  BC  (NBC)  as  the  undirected 


network  having  nodes  E, , ...,  E , N. , N ; for  every  ( j , k)  £ ID, 

I n 1 EQ  D 

there  is  an  arc  (N . , N,)  of  length  b..  between  nodes  N.  and  N,  ; for 
J K J X J K 

every  (i,  j)  £ Ir,  there  is  an  aic  (N.,  E.)  of  length  c..  between  nodes 

u 1 j 1J 

N.  and  E..  We  further  assume  that  the  sets  I-  and  Ir  are  such  that 

1 J D L 

NBC  is  connected,  as  otherwise  DC  decomposes  into  independent  sets  of 


constraints  which  may  be  analyzed  separately. 


Given  a node-path  between  any  two  nodes  f and  f in  NBC,  we  denote 

p q 

the  path  by  P(fp,  and  denote  the  length  of  the  path  by  LP C f ^ , f ) . 

We  define ^(fp,  f^)  to  be  the  length  of  any  shortest  path  in  NBC  between 
nodes  fp  and  f^.  Subsequently,  unless  we  specify  otherwise,  it  should 
be  understood  that  any  path  we  refer  to  is  a simple  path  between  some 
two  nodes  £ and  E . 

P q 

The  following  result  is  established  in  (8j: 

Theorem  2. 1 DC  is  consistent  if  and  only  if 


d(ap,  a^)  <<^(Ep,  E^) , 1 < p < q < n. 


(2.2) 


The  inequalities  (2.2)  are  termed  the  Separation  Conditions  [8], 
since  each  term  on  the  right  specifies  an  upper  bound  on  how  separate 
two  existing  facility  locations  can  be.  Except  when  stated  otherwise, 
we  assume  throughout  this  paper  that  the  separation  conditions  hold,  and 
thus  (equivalently)  DC  is  consistent. 

We  call  a path  P(Ep,  E^)  between  Ep  and  E^  in  NBC  a tight  path  if 
LP(E  , E ) = d(a  , a ).  We  note  that  since  we  assume  DC  is  consistent, 

p q p q 

it  necessarily  follows  if  P(Ep,  E^)  is  a tight  path,  that  LP(Ep,  E^)  =c^( E^ 
Any  path  P(Ep,  E^)  for  which  LP(Ep,  E^)  > d(ap,  a^)  is  called  a 
slack  path,  or  loose  path. 

We  say  that  new  facility  i is  in  a tight  path  if  there  exists  at 
least  one  tight  path  containing  N^.  Every  path  containing  N^  is  slack 
if  there  is  no  tight  path  which  contains  N^. 

The  motivation  for  the  above  terminology  is  due  to  a string  network 
representation  of  NBC.  This  string  network  is  also  useful  for  obtaining 
problem  insights.  When  knots  representing  nodes  Ep  and  E^  are  pulled  as 
for  apart  as  possible,  the  distance  between  the  two  knots  isr^(Ep,  E^). 

If  then  the  string  network  is  placed  upon  the  tree  T, 


V- 


5 


i.e.  the  strings  only  lie  on  arcs  of  T,  a path  is  tight  when  it  is 
necessary  to  pull  the  string  network  tight  in  order  to  place  the  knots 


representing  and  E^  on  a^  and  a^  respectively,  while  a path  is  slack 
if  the  string  path  must  literally  be  slack  when  the  two  knots  are  placed 
to  coincide  with  a and  a . 


p q 

A priori,  one  might  think  that  the  occurrence  of  a tight  path  would 
be  rare.  However,  we  shall  see  in  Sections  4 and  5 that  tight  paths 
occur  in  a quite  natural  way  when  the  separation  conditions  are  used  in 
the  analysis  of  minimax  and  of  efficient  location  problems.  Further, 
the  notion  of  tight  paths  permits  the  specification  of  necessary  and 
sufficient  conditions  for  DC  to  have  a unique  solution. 

Section  3 Preliminary  Analysis  and  Properties  of  Tight  Paths 

In  this  section,  we  establish  some  preliminary  lemmas,  and  then 
present  some  properties  of  tight  paths.  We  will  find  the  following 
three  lemmas  of  use  in  our  analysis.  The  proof  of  Lemma  3.1  can  be 
found  in  [ 7 ] . 

Lemma  3 ■ 1 Given  a,  b,  x £ T with  d(a,  b)  = or  + 3,  d(x,  a)  = a and 
d(x,  b)  = 3,  if  d (y , a)  < a and  d(y,  b)  < f3,  then  x = y. 

Lemma  3 ■ 2 Given  a,  b e T,  d(a,  b)  = or  + f3,  the  inequalities  d(x,  a)  < a, 
and  d(x,  b)  < f3,  are  consistent  if  and  only  if  they  have  a unique 
solution  and  the  inequalities  hold  as  equalities. 

Proof : By  hypothesis  and  the  triangle  inequality,  there  exists  x e T 
where  d(a,  b)  < d(a,  x)  + d(x,  b)  < a + p = d(a,  b),  which  with  the 
hypotheses,  clearly  implies  d(a,  x)  = a and  d(x,  b)  = (3.  Lemma  3.1 
implies  uniqueness.  The  proof  of  the  converse  is  trivial,  and  we  omit 
it. 

Lemma  3.3  Given  x„,  x ..  £ T and  d(xrt,  x . . ) = b.  + b«  + ...  + b , 
0’  r+1  0 r+1  1 2 r+1’ 

the  inequalities 
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I < 1 


are  consistent  if  and  only  if  they  have  a unique  solution  and 
Further, 


d(xi>  W = bi+i  * 1 * 0 r- 


(3.2) 


, x^+^)  i i s •••»  r<  (3.3) 

Proof:  Given  (3.1)  is  consistent,  and  using  the  triangle  inequality  and 
(3.1)  for  any  j , j = 1,  . . . , r, 

d(V  X.)  < a(X„,  X,)  ♦ ...*  dlx.,,  X.)  < b,  x ...  x b.  = 0.  , 


and 


d(xj  ’ W • d(xj>  3tj+l)  + •••+  d(xr  ’ Vl]  ^ bj+i  + •••+  b 


r+1 


Further,  it  is  clear  that  d(x„,  x_.)  = or.  + 8.. 

0 r+1  j 


(3.4) 


Define  vQ  = xQ  , y^  = x^ , and  suppose 

d(7i’  7i+l}  -bi+l  ’ 1 = l»  •**»  c* 

Using  (3.4)  and  the  triangle  inequality,  we  conclude  for  any  j, 

j =1,  ...,  r,  that  d(yQ,  y^)  < cr , and  d(y^ , y^)  < 0 . 

Since  xQ  = yQ  and  x^  = y^ , it  follows  that  d(xQ,  y.)  < or  and 

d(yjt  x^)  £ Pj  • Lemma  3.2  now  implies  y^  = x^  , j = 1,  ...,r.  Repeated 


use  of  the  triangle  inequality  gives  b,  + b-  + ...  +b.,  = d(xn,  x .,) 

I L ITT  l U TiT  1 


< d(xQ,  x1)  + ...  + d(xf,  xr+1)  < bj  + ...  + br+1,  from  which  it 
follows,  using  also  (3.1),  that  (3.2)  is  true.  We  omit  the  (trivial) 
proof  of  the  converse. 

To  show  (3.3),  for  any  i,  note,  in  either  case,  that 


bl  + b2  + •••  + br+l  = d(x0’  Xr+1} 


< d(xQ,  xiM)  + d(x..1,  x.+1)  + d(x.+1,  Xr+1) 


(3.5) 
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But  then 


d(x0’  xi-l)  = bl  + • • • + bi-l* 


d(xi-l*  Vl}  ^ bi  + bi+l, 


d(xi+l>  W = bi+2  + + br+l' 


(3.6) 


Using  (3.5)  and  (3.6),  it  follows  that  d(x.  , , x...)  = b.  + b..,,  which 
along  with  dU^,  x.)  = bi  and  d(x.,  xi+1)  = bi+1,  implies  (3.3). 

We  now  relate  unique  locations  to  tight  paths.  By  definition,  new 
facility  i is  uniquely  located  if  it  has  the  same  location  in  every  fea- 
sible solution  to  DC.  Since  we  later  refer  to  a collection  of  facilities, 
which  contains  possibly  both  existing  and  new  facilities,  being  uniquely 
located,  we  note  that  existing  facilities  are  uniquely  located  by  definition. 

We  now  state  Property  3.1,  the  proof  of  which  follows  from  Lemmas  3.4 
and  3.5  below. 

Property  3.1:  New  facility  k is  uniquely  located  if  and  only  if  N^  lies 
on  at  least  one  tight  path  P(Ep,  E^). 

As  an  immediate  consequence  of  Property  3.1  we  have 
Corollary  3.1:  DC  has  a unique  solution  if  and  only  if  N^  lies  on  at  least 
one  tight  path  in  NBC  for  k = 1,  . . . , m. 

We  now  present  Lemmas  3.4  and  3.5. 

Lemma  3.4:  If  N^  lies  on  at  least  one  tight  path  P(Ep,  E^)  in  NBC,  then 
new  facility  k is  uniquely  located. 

Proof : Index  the  nodes  in  the  tight  path  P(Ep,  E^)  as 

(fQ,  fj,  f2,  ...,fr,  f^j),  where  fQ  = Ep,  f^  = Eq  and  each  f ^ , 2 < j < r, 
is  an  N or  E node.  (We  note  that  since  no  two  E nodes  are  adjacent  in 
NBC,  no  two  adjacent  nodes  in  (£q,  f^,  ...,  f are  E nodes.) 
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Letting  b.  be  the  length  of  arc  (f.,,  f.),  j = 1,  ....  r+ 

J J ^ J 

since  P(E  , E ) is  a tight  path  we  have 

p q 

d(V  aq)  = bj  + b2  + ...  + br+1. 

Since  DC  is  consistent,  there  exists  a location  vector 

X = (Xj,  x^)  which  satisfies  the  distance  constraints. 

Let  j be  arbitrary,  0 < j < r. 

If  fj  = N(j)  aQd  fj+i  = E(j+i)’  then  e JC  aad 

d(x(j)’  a(j+D}  - C(j),(j+1)  ' bj+l‘ 

If  fj  = E(j)  and  fj+l  = N(j+1)’  then  + e Jc  and 

d(aCj) ’ xU+i))  = d(xCj-*-D ’ a(j)}  - C(j+D,(j)  5 bj+r 

If  fj  = N(j)  aad  fj+l  = N(j+1)’  thea  e IB  and 

d(x(j)’  X(j+D)  - b(j),(j+n  5 bj+r 

But  then,  from  (3.7),  (3.8),  (3.9),  and  (3.10),  the  hypotheses  of 
Lemma  3.3  are  satisfied  and  so  new  facility  k is  uniquely  located. 

The  conclusion  now  follows. 

We  next  establish  the  converse  of  Lemma  3.4. 

Lemma  3.5:  If  new  facility  k is  uniquely  located  then  lies  on  3t 
least  one  tight  path  P(Ep,  Eq)  in  NBC. 

Proof:  We  shall  show  equivalently  that  if  every  path  containing  N^  is 
slack,  then  there  exist  solutions  to  DC  with  distinct  locations  of  new 
facility  k,  where,  without  loss  of  generality,  we  let  k = 1.  Assume  all 
paths  through  Nj  are  slack.  It  is  then  the  case  that  any  such  path  with 
"least"  slack  has  positive  slack.  That  is.  there  are  existing  facility 
nodes,  taken  to  be  E^  and  E2  without  loss  of  generality,  and  a path 
P12  = P(E1,  Ej)  which  contains  Nj,  such  that,  with  dj,  = dUj,  a?), 

8 = EP12  - dl2  > 0 and,  for  any  E and  E and  any  path  £ = $(E  , E ) 

containing  Nj,  with  djt  s d(as,  *t),  it  is  true  that 


i; 

(3.7) 

(3.8) 

(3.9) 

(3.10) 
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L?  > 6 + d 
st  - p st 


(3.11) 


Since  has  least  slack,  it  is  a shortest  path  through  Nj 

joining  and  and  thus  LP^  ^^(Ej,  Nj)  +<5^(E2,  N^).  As  a notational 

simplification,  we  let  JL  =<^)(Ei,  N^,  i = 1,  2. 

Let  DC  , i = 1,  2,  be  the  set  of  distance  constraints  consisting 

of  DC  in  addition  to  the  constraint  d(x^,  a^)  < y^.  We  note  that  any 

feasible  solution  to  DC,  i = 1,  2,  is  a feasible  solution  to  DC.  In 

what  follows  we  show  that  by  proper  choice  of  y^,  DC^  is  consistent, 

i = 1,  2,  and  that  if  X1  is  feasible  to  DC^  and  is  feasible  to  DC,, 

2 2 

then  x1  and  x^  are  distinct,  where  x^  is  the  first  component  of  X1, 
i = 1,  2. 


Let  NBC.  be  the  network  associated  with  DC..  We  note  that  NBC 
i ii 

differs  from  NBC  only  in  that  it  contains  the  additional  arc  between 
nodes  E.  and  N,  of  length  y..  Label  this  additional  arc  a.. 

ill  i 

Let  E and  E be  any  two  existing  facility  nodes  in  NBC.  and 

S L 1 

let  P be  a shortest  path  in  NBC.  of  length  LP  „ , joining  E and  E . 

S L 1 St  SC 

We  shall  show  that  if  y^  is  positive  and  if 

- P , (3.12) 

then, 


LP  > d , 
st  - st’ 

and  thus,  due  to  Theorem  2.1,  DC.  is  consistent.  We  shall 

l 


(3.13) 


show  that  (3.13)  holds  when  i = 1.  Due  to  symmetry,  the  case  when  i = 2 
will  follow. 


If  Pst  contains  only  arcs  which  are  arcs  of  NBC,  i.e.,  P does  not 
contain  a ^ , then  (3.13)  is  certainly  true  since  the  separation  conditions 
hold  for  DC.  Thus  suppose  Pgt  contains  a (and  hence  contains  N^ ) . P 
can  be  decomposed  as  Pjt  = (P^,  or^,  P^),  where  Pj  and  P,  are  (possibly 
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f 


empty)  paths  joining  E and  E1 , and  N.  and  E , respectively.  Since  P 
is  a shortest  path,  neither  P1  nor  P^  contains  a^,  and  thus  P^  and  P^ 

(when  non-empty)  contain  only  arcs  which  are  arcs  of  NBC.  The  length  of 
Pgt  is  = LPj  + Yj  + LP2,  and  using  (3.12)  we  obtain 

!Pst  £ LPi  +Z±  - 3 + lp2’  or 

LPst  + 0 > LP1  + £j  + LP2.  (3.14) 

Noting  that  the  right  hand  side  of  (3.14)  is  the  length,  1$  , of  a 

path,  ? , in  NBC  between  E and  E which  contains  N, , (3.13)  follows 

St  s c 1 

from  (3.11)  and  (3.14).  (We  remark,  to  motivate  the  requirement  (3.12), 
that  if,  say,  Yj  < 2^  - 3,  then  + 2^  < 2^  + 2^  - $ = di2  : but 
Yj  + 22  is  t^e  length  of  a path  in  NBC^,  not  entirely  contained  in  NBC, 
and  so  Theorem  2.1  would  imply  DC^  is  inconsistent. ) 

Define  the  positive  quantity  e = min{2^,  2^,  d32},  an<*  ^et  ^i  be 
the  positive  quantity,  which  we  note  satisfies  (3.12),  defined  as 
Y^  = max{£/3,  JL  - {3} , i = 1,  2.  Since  DC^  is  consistent,  let  X1  be  a 
feasible  solution  to  DC^,  and  let  x*  be  the  first  component  of 
X1,  i =1,  2. 

For  x^  we  have  d(a1?  Xj)  < Yj  = max{  / ^ * * P}»  and  so  the 

triangle  inequality  gives  d(a2,  xj)  > d^  - d(a1,  x*)  > d12  - 
max{£/3,  - 3}  = min{dj2  - S/3,  dJ2  - + 0}.  Noting  that 

d^2  “ 2^  + 0 = 2^  we  obtain 

d(a2,  xj)  > min{d12  - e/3>  2^} . (3.15) 

2 

For  x^  we  have 

d(a2,  x^)  < y2  = «ax{£/3,  22  - 0} . (3.16) 

Using  the  definition  of  e,  it  follows  that  the  right  hand  side  of  (3.15)  is 

1 2 

strictly  larger  than  the  right  hand  side  of  (3.16)  and  thus  x1  and  x^  are 
distinct.  Noting  that  xt  i = 1,  2,  is  feasible  to  DC  completes  the  proof. 


1 


11 


We  can  now  draw  some  additional  conclusions  about  tight  paths. 
Property  3.2:  If  P(Ep,  E^)  is  a tight  path,  then  the  nodes  representing 
facilities  in  the  path  occur  with  the  same  ordering  and  spacing  in  the 
path  as  do  the  locations  representing  the  facilities  in  L(ap,  a^) . 
Further,  every  facility  represented  by  a node  in  P(E  , E^)  is  uniquely 
located. 


Proof : The  proof  is  an  immediate  consequence  of  Lemmas  3.3  and  3.4. 

We  emphasize  the  fact  that  the  knowledge  of  a tight  path  immediately 
identifies  the  location  of  every  new  facility  represented  by  a node  on 
the  tight  path.  For  example,  if  P(E1>  E-j)  = (E^  N^  N2,  E2,  Nj,  E3)  is 
a tight  path,  with  cn  = 2,  b12  = 3,  = 1,  c32  = 4,  c33  = 2,  then 

d(a^,  a3)  = 12,  and  so  new  facility  1 would  be  located  in  L(a^,  a3)  such 
that  d(x, , a^)  = 2,  new  facility  2 would  be  located  in  L(a1,  a3)  such 
that  d(x2,  a^  = 5,  and  new  facility  3 would  be  located  in  L^,  a3)  such 
that  d(x3,  a1)  = 10. 

We  now  consider  the  problem  of  determining  when  an  arc  lies  on  a 
tight  path.  As  an  arc  lies  on  a tight  path  if  and  only  if  it  is  not  the 
case  that  all  paths  containing  the  arc  are  slack,  we  consider  the 
equivalent  problem  of  determining  when  an  arc  lies  only  on  slack  paths. 
Property  3.3:  Let  DC  be  consistent.  Let  ( f ^ , f ) be  any  arc  in  NBC,  of 
length  e_ , whose  length  is  reduced  by  some  positive  amount  £.  Let 
DC£(NBC£)  be  the  distance  constraints  (network)  obtained  from  DC (NBC)  by 
replacing  e^.  by  e^  - £. 

(a)  Every  path  containing  (f^,  f ) in  NBC  is  slack  if  and  only  if  £ can 
be  chosen  (with  £ > 0)  so  that  DC  is  consistent. 

(b)  Whenever  every  path  containing  ( f ^ , f ^ ) is  slack,  £ can  be  chosen 

(with  £ > 0)  so  that  DC  is  consistent  and  at  least  one  of  the  following 

& 

is  true: 
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(i)  at  least  one  path  in  NBC  containing  (f . , f.)  is  tight; 

t 1 J 

(ii)  the  length  of  ( f . , f.)  in  NBC  can  be  reduced  to  zero. 

1 J C 

Proof:  Assume  e can  be  chosen  (with  £ > 0)  so  that  DC  is  consistent. 

£ 

Theorem  2.1  then  implies  the  separation  conditions  hold  for  DC  . If  at 

least  one  path,  say  P(Ep,  E^),  in  NBC  which  contains  (f^,  f ^ ) is  tight, 

then  d(ap,  a^)  = LP(Ep,  E^) , so  that  reducing  e^  by  £ would  give 

EP£(Ep,  E^)  = LP(Ep,  E^)  - £ < dCa^,  a^).  But  LP£(Ep,  E^)  is  the  length  of 

P(E  , E ) in  NBC  , and  hence  a separation  condition  for  DC  is  violated 
p q t £ 

for  any  £,  £ > 0,  giving  a contradiction.  Hence  every  path  containing 

(f.,  f.)  in  NBC  is  slack, 
i*  J 

Conversely,  assume  every  path  containing  (f^  f ^ ) in  NBC  is  slack, 


so  that 


LP(E  , E ) - d(a  , a ) > 0 
P q P q 


(3.17) 


for  every  path  P(Ep,  E^)  containing  (f^,  f ^ ) , and  define  £'  to  be  the 
minimum  of  the  left  side  of  (3.17)  over  all  such  paths,  giving  £'  >0. 
Let  £ = min(e' , e„)  > 0.  Since  any  path  P(Ep,  E^)  containing  ( f ^ , f^) 
has  its  length  reduced  to  LPe.(Ep,  E^)  = LP(Ep,  E^)  - £,  the  inequality 


d<V  V - ‘VV  V 


£ < L-‘(E  , E ) - d(a  , a ) . 
- p’  q P’  q' 


(3.18) 

H «•  if  4 

is  equivalent  to 

(3.19) 

The  definition  of  £'  and  £ < e’  imply  (3.19)  is  true,  and  so  (3.18)  is 

true.  Further,  if  e = S',  then  (3.19)  holds  as  an  equality  for  some 

P(Ep,  E^) , implying  its  length  in  NBC£  is  zero,  giving  (b)  - (i).  The 

case  £ = e^  gives  (b)  - (ii).  (Parenthetically,  we  note  that  the 

largest  £ can  be  and  still  have  DC£  consistent  is  thus  clearly 

max(£',  e..).)  As  any  path  P(E  , E ) in  NBC  not  containing  (f.,  f.) 

lj  P 4 E 1 J 

has  its  length  unchanged  from  that  in  NBC,  DC  consistent  and 

Theorem  2.1  implies  the  path  satisfies  its  separation  condition  for  DC£. 


13 


Hence,  with  (3.18),  the  separation  conditions  for  DC£  are  satisfied,  so 
Theorem  2.1  implies  DC£  is  consistent,  completing  the  proof. 

Section  4:  A Minimax  Location  Problem  Given  m new  facilities  at  locations 
Xj,  ...,  to  be  determined,  existing  facilities  at  known  locations  a^, 
...,  aQ,  and  index  sets  Ig,  1^  as  before,  the  minimax  location  problem 
may  be  stated  as  follows: 

(PMM)  minimize  z 
subject  to 

vjk  d(Xj , xfc)  < z , (j , k)  £ Ig 

wtj  d(xi , a^)  < z , (i,  j)  e Ic. 

Here  the  v.,  and  w are  given  positive  weights,  and  the  problem  becomes 

jx  1J 

one  of  locating  the  new  facilities  so  as  to  minimize  the  maximum  of  the 

weighted  distances.  To  define  the  network  BC  of  interest  for  this 

problem,  define  b.,  = 1/v.,  for  ( j , k)  e I_,  and  define  c. . = 1/w. . for 
jk  jk  a ij  ij 

(i,  j)  £ Ic. 

The  following  result  for  this  problem  is  proven  in  [8]: 

Lemma  4.1:  The  minimum  objective  function  value  for  (PMM)  is  given  by 
z*  = max[d(ap,  aq)/^(Ep,  Eq)  : 1 < p < q < m] . 

As  is  pointed  out  in  [8],  one  way  to  solve  the  minimax  problem  is  to 
use  Lemma  4.1,  set  z = z*  in  the  constraints  of  the  problem,  and  use  the 
SLP  algorithm  to  construct  a feasible  solution  to  the  constraints;  any 
such  feasible  solution  is  optimal. 

Our  interest  here  is  to  explore  relationships  between  Lemma  4.1 
and  tight  paths.  To  this  end,  it  will  be  convenient  to  denote  by  NBC(z) 
the  network  obtained  from  NBC  by  multiplying  every  arc  length  in  NBC  by 
z.  We  now  state 
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Property  4.1;  Let  (X,  z)  be  a feasible  solution  to  (PMM). 

(a)  (X,  z)  is  an  optima  feasible  solution  to  (PMM)  if  and  only  if  at 
least  one  path  in  NBC(z)  is  tight,  that  is,  for  soae  P(Ej , E^) . 

d(ajt  ak)  = z LP(E.,  Efc). 

(b)  For  any  such  tight  path,  the  facilities  whose  nodes  lie  on  the  path  are 
uniquely  located,  and  their  locations  have  the  same  ordering  and  spacing  in  T 
as  their  nodes  have  in  the  corresponding  path  in  NBC. 


Proof : (a)  Suppose  P(E.,  E.)  is  a tight  path,  so  that 

J * 


d(..,  .k)>zLP(E.,Ek).  ( 

The  definition  of^(Ej,  E^)  and  z > 0 give 

z LP(Ejf  Ek)  > z ^(Ej,  Ek),  ( 

while  z*  the  ainimuai  objective  function  value  (by  Learn  4.1)  gives 
z ^(Ejf  Ek)  > **  ^(E j , Ek).  ( 

Hence  LeasM  4.1  gives 

z*^(E.,  Ek)  > ld(a.,  ak)/^(Er  Ek)]<^(E.,  Efc) 


(4.1) 


(4.2) 


(4.3) 


= d(aj , ak).  (4.4) 

Thus  we  conclude  the  inequalities  (4.1)  through  (4.4)  all  hold  as 
equalities,  implying  z = z*,  and  hence  (X,  z)  is  an  optiaua  feasible 
solution. 

Conversely,  suppose  (X,  z)  is  an  optiaua  feasible  solution.  Leaai 
4.1  then  gives,  for  soae  P(E^,  Ek), 

z = z*  = aax[d(ap,  *q)/c^(Ep,  Eq)  : 1 < p < q < n] 

* "'•j-  V 

= d(Sj , •k)/LP(Ej , Ek), 
iaplying  P(E^,  Ek)  is  a tight  path  in  NBC(z). 

(b)  The  proof  is  iaaediate  froa  (a)  and  Properties  3.1  and  3.2. 
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We  now  present  an  example.  Suppose  a tree  T is  given  as  shown  in 
Fig.  1,  with  three  new  facilities  and  three  existing  facilities 
(at  ap  a2,  and  83).  Ig  = {(1,  2),  (1,  3),  (2,  3)},  Ic  = {(1,  1), 

(2,  2),  (3,  3)},  v12  = 1,  v13  = 10,  v23  = 1,  wn  = 10,  w22  = 2,  and 

w,,  = 5.  Fig.  2 illustrates  NBC  with  b..=  1/v..  and  c. . = 1/w... 

33  0 jk  jk  ij  ij 

It  is  readily  verified  that  ^(£3,  E2^  = E3)  = *4, 

^(E2,  E3)  = 1.7,  d(a^ , a2)  = 5,  d(a1>  a3)  = 6,  and  d(a2,  83)  = 7. 

Thus  Lemma  4.1  gives  2*  = raax[5/.1.6,  6/. 4,  7/1.7]  = 15  * d(aj,  ajJ/^CEj,  E3). 
Hence,  in  NBC(z*),  Property  4.1  implies  the  path  (E^,  N^ , N3,  E3)  is  tight, 
with  total  length  ( . 1 ) ( 15 ) + (.1)05)  + (.2)05)  = 6 = d(a^,  aj),  and  so 
new  facility  1 and  new  facility  3 are  uniquely  located  at  x*,  x|  respectively 
in  L(a^,  83),  where  x*  is  the  point  in  Ka^,  83)  such  that  d(a^,  x*)  = 1.5, 
and  x|  is  the  point  in  L(a^,  a 3)  such  that  d(a^,  x^)  =3.  It  is  easy  to 
verify  that  N2  lies  only  on  loose  paths  in  NBC(z*),  and  hence  new  facility 
2 is  not  uniquely  located.  In  fact,  the  inequalities  2d(x|,  a2)  < 15, 
ld(x*,  x*)  < 15,  and  ld(x*,  x*)  < 15  permit  x*  to  be  any  point  in  the 
tree  T.  This  example  illustrates  the  fact  that  the  knowledge  of  tight 
paths  in  NBC(z)  permits  one  to  determine  immediately  some  new  facility 
locations  for  (PMM) , as  well  as  identify  those  facilities  which  have 
"critical”  locations,  in  the  sense  that  a change  in  their  locations 
would  cause  z*  to  change  as  well. 
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Section  5 Efficient  Location  Vectors 


As  before  we  suppose  existing  facility  locations  a^,  aQ  to  be 

given,  and  that  new  facility  locations  Xj , . ..,  x^  are  of  interest. 

Given  any  vector  of  new  facility  locations  Z = (z^  zm),  we  denote 

by  D(Z)  the  vector  having  the  entries  d(z.,  z.  ) for  ( j , k)  e In,  and 

j k a 

d(z.,  a.)  for  (i,  j)  e Ir,  where  ID  and  I-  are  given  pairs  of  facility 
1 J v-  o L 

indices.  A location  vector  Z is  dominated  if  there  exists  a location 

vector  Y such  that  D(Y)  < D(Z)  and  D(Y)  t D(Z).  A location  vector  Z 

which  is  not  dominated  is  called  efficient.  Thus  Z is  efficient  if  and 

only  if  D(Y)  < D(Z)  implies  D(Y)  = D(Z).  Hence,  whenever  Z is  efficient, 

if  a location  vector  X is  given  such  that  some  entries  in  D(X)  are 

(strictly)  less  than  the  corresponding  entries  in  D(Z),  then  it  must  be 

true  that  at  least  one  entry  in  D(X)  is  (strictly)  greater  than  the 

corresponding  entry  in  D(Z).  One  can  consider  the  problem  of  finding 

efficient  location  vectors  as  a multiobjective  optimization  problem, 

with  one  optimizer  for  every  entry  in  D(Z);  the  optimizers  can  agree 

that  dominated  location  vectors  are  not  of  interest,  thus  leaving  the 

efficient  location  vectors  to  consider.  Even  if  there  is  only  a single 

optimizer,  and  his  objective  function  is  strictly  increasing  in  each 

entry  of  D(Z),  then  every  location  vector  which  minimizes  his  objective 

function  is  clearly  efficient,  so  that  knowledge  of  efficient  location 

vectors  may  facilitate  a sensitivity  analysis. 

As  an  example  of  such  an  objective  function,  suppose  positive 

''weights”  v and  w.  . are  given,  and  define  f(Z)  by  f(Z)  = 

J * • J 

2{yjk  d(zj»  zk)  : (J.  k)  6 IB)  + SNij  d(zi,  *j)  : (i,  j)  £ ICJ- 

Picard  and  Ratliff  [17]  have  recently  presented  a means  of  finding  new 
facility  locations  in  a tree  to  minimize  f(Z),  and  discuss  related 
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literature,  which  is  extensive.  Whenever  a single  optimizer  has  a 
continuous  objective  function  (such  as  f(Z)),  then  the  compactness  of 
T and  the  extreme  value  theorem  asserts  there  will  always  exist  a loca- 
tion vector  optimizing  his  objective  function,  in  turn  implying  that 
efficient  location  vectors  always  exist  for  nonpathological  problems. 

A further  discussion  of  multiobjective  location  problems  on  networks  can 
be  found  in  [11],  [12],  [13],  and  [16]. 

Given  a location  vector  Z,  we  let  U = D(Z)  and  define  the  distance 
constraints  of  interest  by  D(X)  < U,  where  the  entries  in  U define  the 
bjk  and  c^  by  = d(z^,  zk)  for  ( j , k)  s Ig,  and  c^  = d(zi,  a^)  for 
(i,  j)  £ Ij,.  We  use  the  b^k  and  c^.  to  define  NBC  in  the  customary 
manner.  As  before,  we  may  assume  NBC  is  connected,  for  otherwise  the 
problem  of  finding  efficient  location  vectors  decomposes  into  independent 
subproblems.  Further,  we  note  that  DC  is  always  consistent,  as  Z is 
certainly  feasible  to  DC,  and  hence,  by  Theorem  2.1,  the  separation 
conditions  are  always  satisfied.  For  convenience,  for  any  location 
vector  Z,  we  denote  by  A^(Z)  the  collection  of  locations  of  uniquely 
located  facilities  whose  nodes  are  adjacent  to  N^  in  NBC.  We  denote 
by  H(A^(Z)]  the  convex  hull  of  A^CZ),  the  imbedding  of  the  smallest 
subtree  of  T spanning  all  the  elements  of  A^(Z). 

With  the  above  definitions  we  can  present  a family  of  equivalent 
conditions  for  a location  vector  Z to  be  efficient. 

Property  5.1:  Given  a location  vector  Z used  to  define  DC  and  NBC,  the 
following  are  equivalent: 

(a)  Z is  efficient, 

(b)  Each  is  in  at  least  one  tight  path  in  NBC, 

(c)  Z is  the  unique  solution  to  DC, 

(d)  z.  e H[A.(Z)]  for  i = 1 m. 
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Proof:  The  equivalence  of  (b)  and  (c)  is  a direct  consequence  of 
Property  3.1  and  the  fact  that  Z is  always  a feasible  solution  to  DC, 
while  (c)  clearly  implies  (a).  To  show  (a)  implies  (b),  suppose  some 
is  not  in  a tight  path.  Property  3.3  then  implies  some  entry  in  U = D(Z) 
can  be  reduced  and  the  resultant  distance  constraints  will  still  have  a 
feasible  solution,  say  Y.  But  then  clearly  D(Y)  < D(Z)  and  D(Y)  ^ D(Z) , 
contradicting  the  fact  that  Z is  efficient.  Hence  (a),  (b),  and  (c)  are 
equivalent.  It  can  be  seen  that  the  proof  will  be  complete  if  we  show 
(b)  implies  (d)  and  (d)  implies  (c) . 

To  show  (b)  implies  (d) , suppose  is  in  some  tight  path  P.  Let 
f^  and  f2  be  the  nodes  adjacent  to  JT  in  P,  so  that  ((f  , }T) , (N^,  f^)) 
is  a subpath  of  P.  Since  f^  and  f2  are  in  the  tight  path  P,  by  Property 
3.2  the  facilities  represented  by  f^  and  f^  are  uniquely  located.  We 
may  let  y^  and  y2  denote  the  unique  locations  of  f^  and  f^  respectively. 

Thus  it  is  clear  that  y^  and  y2  are  elements  of  A^(Z).  By  Property  3.2, 
zi  e LCyj,  y2),  and  by  definition  of  the  convex  hull,  L (y ^ , y2)  C.  HfA^CZ) ) . 
Thus  it  follows  that  z ^ e HCA^CZ)).  To  show  (d)  implies  (c),  suppose 
z^  e H[A^(Z)]  and  let  f^  and  f2  be  nodes  adjacent  to  N^  in  NBC,  where 
f^  and  f2  represent  facilities  with  unique  locations  y^  and  y2>  respectively, 
such  that  zi  £ L(ylt  y2)<CT.  Thus  d(yx,  y2)  = d(y1?  z /)  + d(zi>  y2) . 

Now  for  any  feasible  solution  X to  DC  we  know  dCy^  x^)  < d(y1,  z^  and 
d(y2,  x.)  < d(y2,  z.).  But  then  because  f^  and  f 2 are  uniquely  located, 

Lemma  3.2  implies  x^  = z^,  for  i = 1 , . . . , m.  Hence  X = Z,  so  Z is  the 
unique  solution  to  DC,  completing  the  proof. 

As  an  example,  suppose  we  again  have  three  existing  and  three  new 
facilities,  and  use  the  tree  of  Fig.  1.  Suppose  the  entries  of  D(Z)  are 
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given  by  d(zr  a^,  d(zr  a3),  d(z2,  a1),  d(z2,  a3),  d(z3>  a^,  d(z3>  a2) 
d(z3>  z2) , d(z1,  z3),  and  d(z2>  z3).  Suppose  the  location  vector  Z is 
as  shown  in  Fig.  3.  Fig.  3 and  the  above  distances  gives  NBC  as  shown 
in  Fig.  4.  As  the  path  (E^,  N^  E3)  in  NBC  has  length  6 = dCa^  a3), 

Nj  is  on  a tight  path.  Further,  the  path  (Ej,  N2>  E3)  also  has  length 
6 = d(a^,  a3),  so  N2  is  on  a tight  path.  However,  it  can  be  verified 
that  every  path  containing  N3  is  slack,  so  the  vector  Z = (z^,  z2,  z3), 
is  dominated.  It  is  easily  verified  however  that  if  z3  is  changed  so 
that  either  z 3 £ L(a^,  a2)  or  z3  £ LCz^,  z2) , that  the  new  location  vector 
will  be  efficient,  as  N3  is  on  at  least  one  tight  path,  or,  equivalently, 
z3  £ H({alf  a2,  z2,  z2}). 

As  a final  application  of  the  idea  of  efficiency,  let  us  choose  an 

efficient  location  vector  problem  and  an  efficient  location  vector  Z for 

which  U = D(Z)  will  have  m(m-l)/2  + mn  entries.  Suppose  we  subtract  any 

positive  quantity  £ from  an  arbitrary  entry  in  U to  obtain  a vector  U£, 

giving  U£  < U and  U£  t U.  The  distance  constraints  D(X)  < U£  must  then 

be  inconsistent,  as  else  X would  dominate  Z.  Now  given  the  constraints 

D(X)  < U£,  if  we  do  not  know  which  entry  in  U has  been  reduced  by  £ to 

obtain  U , then  any  algorithm  we  apply  to  check  if  D(X)  < U is  consistent 

must  examine  all  of  the  m(m-l)/2  ♦ mn  constraints,  since  any  one  of  the 

constraints  can  cause  inconsistency.  Hence  any  algorithm  to  determine 

if  D(X)  < U£  is  consistent  is  at  least  of  order  m + mn.  Thus  we  have 

Property  5.2:  Any  algorithm  to  determine  whether  or  not  the  distance 

constraints  of  Section  2 are  consistent  is  (for  the  worst  case)  at  least 
2 

of  order  m + mn. 

As  is  pointed  out  in  [8],  the  "Sequential  Location  Procedure"  presented 

2 

in  [8]  to  check  the  consistency  of  DC  is  of  order  m + mn,  and  thus 
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Property  5.2  implies  it  is  in  fact  an  algorithm  of  lowest  order  for 
checking  consistency. 
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